Exploring a deep meridional flow hypothesis for a circulation 
dominated solar dynamo model 

G. A. Guerrero* J. D. Munozf Elisabete M. de Gouveia Dal Pino* 

February 2, 2008 

Abstract 

Circulation-dominated solar dynamo models, which employ a helioseismic rotation profile 
and a fixed meridional flow, give a good approximation to the large scale solar magnetic phe- 
nomena, such as the 11-year cycle or the so called Hale's law of polarities. Nevertheless, the 
larger amplitude of the radial shear dQ/dr at the high latitudes makes the dynamo to produce 
a strong toroidal magnetic field at high latitudes, in contradicti on with the observati o ns of the 
sunspots (Sporer's Law). A possible solution was proposed bv lNandv fc Choudhuril i2002T) in 
which a deep meridional flow can conduct the magnetic field inside of a stable layer (the ra- 
diative core) and then allow that it erupts just at lower latitudes. Although they obtain good 
results, this hypothesis generates new problems like the mixture of elements in the radiative 
core (that alters the abundance of the elements) and the transfer of angular momentum. We 
have recently explored this h ypothesis in a different approxi mation, using the magnetic buoy- 
ancy mechanism proposed bv lDikpati fc Charbonneaul il999T) and found that a deep meridional 
flow pus hes the maximum of the tor oidal magnetic field towards the solar equator, but, in con- 
trast to lNandv fc Chou dhuril i2002l) . a second zone of maximum fields remains at the poles 
dCuerrero fc Munozl I2004T) . In that work, we have also introduced a bipolytropic density pro- 
file in order to better reproduce the stratification in the radiative zone. We here review these 
results and also discuss a new possible scenari o where the tachocline has an ellipsoidal shape, 
following early helioseismologic observations iCharbonneau et all [l999), and find that the 
modification of the geometry of the tachocline can lead to results which are in good agreement 
with observations and opens the possibility to explore in more detail, through the dynamo 
model, the place where the magnetic field could be really stored. 



1 Introduction 

Over the last few years, large scale solar magnetic field features, such as the 11 years sunspots 
cycle or polarity inversions of the magnetic field, have been successfully explained by solar 
dynamo models that include a solar like differential rotation and meridional circulation as the 
main ingredients. Since helioseismic measurements have indicated the presence of a tachocline, 
it has been a common believe that the dynamo action takes place in this thin layer with 
substantial radial shear. In the tachocline, a poloidal magnetic field is stretched by the solar 
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rotation such that a belt of strong toroidal magnetic field is formed around the solar equator, 
theoretical arguments suggest that this field is not homogeneous but mainly concentred in 
magnetic flux tubes surrounded for less magnetized plasma. The strong magnetic pressure 
can make these flux tubes rise to the surface by the action of magnetic buoyancy. There, 
the poloidal field is possibly regenerated by the decay of the bipolar magnetic regions. This 
poloidal field is then advected by meridional circulation towards the solar poles and then to the 
deep layers, where the full cycle will be completed. It is not completely understood yet what 
is the profile of the meridional circulation, but its crucial role on the dynamo process has been 
recently recognized and the models that include it have been named circulation dominated 
solar dynamo (CDSD) models. 

Most of the models that consider a differential rotation profile derived from helioseismic 
inversions, present an interesting problem, the radial shear dQ./dr (where f2 is the angular 
velocity and r is the radial distance) is larger at higher latitudes than at lower ones so that 
a strong toroidal field is expected to be generated in regions closer to the poles and therefore 
sunspots should appear at higher latitudes, contradicting the observations (Sporer's Law). 
Apparently this was an unavoidable problem and no ch ange in the parameters cou ld solve it 
iDikpati fc Charbonneaulll999l . lKuker et alll200ll) until lNandv fc Choudhuril (I2002D proposed 
a new possible scenario in order to reduce the amplitude of the toroidal field at higher latitudes. 
They suggested that the meridional flow could go a little deeper than in the previous models. 
With a flow going below the tachocline, the toroidal magnetic field is stored in a more stable 
layer and advected by meridional circulation to lower latitudes from where it emerges to the 
tachocline and the con vective zones and undergoes magnetic buoyancy. This model (see also 
IChatteriee et all i2004) ') results a magnetic field distribution which is in qualitative agreement 
with the observations of the sunspots distribution, though it may affect the mixture of the 
elements and the observed abundances, and also the angular momentum transfer, but these 
specific potential problems will not be addressed in the present work. 

Th ere is also another fundamental difference in the way by which the model of lNandv fc Choudhuril 
( 2002() and the previous ones treat the buoyant process. One of the aims in these models is 
to produce a flux tube of strong toroidal field (B c = 10 G) in a subadiabatic layer that 
becomes buoyant unstable and emerges to the surface in a time which is much smaller (a 
month) than the characteristic time of one magnetic cycle (of the order of ten years). While 
iDikpati fc Charbonneaul <ll999Tl used a nonlocal treatment, making the poloidal field to regen- 
erate clos er to the surface proportion ally to the toroidal field at the base of the solar convec- 
tion zone. lNandv fc Choudhuril i2002l) employed a numerical procedure in which whenever the 
toroidal magnetic exceeded B c at the base of the solar convection zon e, a fraction / of it is 
made artificially to eru pt to the surface at time in tervals r = 8.8 x 10 5 iNandv fc Choudhuril 
1200 ll) . In a recent work lGu errer^^^M^noj ||200^), we ha ve developed a CDSD model using a 
solar differential rotation like iNanch^^Chouofliui^ ( 2002 ) , but employed the magnetic buoy- 
ancy mechanism used by IDikpati fc Charb^nneau W9i ). We have a lso considered a more 
realistic bipolytropic density distribution jPin^6nfcCalvo-Mozol . 1200 ll) in order to better re- 
produce the subadiabatic stratification of the radiative zone. We have found that properties 
such as the 11 years cycle, the inversion of polarity, and the equatorward migration of branches 
of strong toroidal field and the poleward migration of a weak poloidal field, are correctly re- 
produced, however even using a deep meridional flow, a strong toroidal field still remains in 
the poles. The introduction of the bipolytropic profile has not modified either the morphology 
of the butterfly diagrams or the branches of strong toroidal field in the poles. 

In the present work, we have tried to improve our model by modifying the geom e try o f 
the tachocline taking into account the early helioseismic results of ICharbonneau et alJ il999t) 
according to which the tachocline may have a prolate shape. In section we present the math- 
ematical formalism of our model. The numerical details can be found in (G uerrero fc Munoa 
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12004) . Section JUcont ains the main results of our CDSD model with a prolate tachocline and, 
finally, in we draw our conclusions. 



2 The Model 

Our code solves the MHD induction equation which governs the evolution of the solar dynamo: 

V x (U x B) + 77V B . (1) 



dt 

By assuming spherical symmetry, the magnetic and velocity fields can be writen as 

B = J?M,i) + Vx (A(r,6,t)), (2) 
U = u(r,8) +r sin 9Q(r,6) , (3) 

where B(r,6,t) and V x (A(r, 8, t)) correspond to the toroidal and the poloidal components 
of the magnetic field, respectively; Q is the angular velocity, u = u r + ug is the velocity in 
the meridional plane, and 77 is the magnetic diffusivity. Replacing equations @ and (J^J in the 
induction equation Q and separating the poloidal and toroidal components of the magnetic 
field, we obtain 

^- + l(u-X7)(sA)= Vp (X7 2 -±)A + S 1 (r,e,t) , (4) 

^ + ;[|:(r«rB) + ^(« fl S)] = (Bp.V)n (5) 

-V-qt x V x B + r^V 2 - , 

where s=r sin 8 and B P =V x A. The values and profiles of r\ v and r\ t will be disscused in the 
section l2~4l In the right side of eq. Si(r, 8, B^) is a source term (see below). 



2.1 Differential rotation 

The solar differential rotation inferred from helioseismology can be described by the following 
equation 

n(r,9) = n c + ±[l + erf(2^±)](n s (6)-n c ) . (6) 

Here, fl s (d)=QEq + ci2 cos 2 8 + as. cos 4 8 is the latitudinal differential rotation at the sur- 
face and erf(x) is an error function that confines the radial shear to a tachocline of thick- 
ness di=0.05-Re. In this expression, a rigid core rotates uniformly with angular velocity 
0:/27r=432.8. The other values are f2 Eq /27r=460.7, a 2 /27r=-62.9, a 4 /27r=-67.13 nHz, and 
r c =0.7R Q . 



2.2 Meridional circulation 

As in iM. Dikpatil 11994 lA. R. Choudhuri Il995l iNandv fc Choudhuri I2002T) . we assume a 
single convection cell for each meridional quadrant 

p(r)u = Vx [1>(r,6)e<,] , (7) 
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where ip is the stream function given by 



Ursine = ( r -i^ 0Sm [ g?_- || ] (8) 

x (l_ e -' 9 ^ fl£ )(l- e ^( e -/ 2 )) 

x e [(r-ro)/r] 2 

and p is the density profile for an adiabatic gaseous sphere with a specific heat ratio coefficient 
7=5/3 (m=f.5). Then 

p(r) = C(^-0.95)" 1 , (9) 

wher e we chose C=3.60 x 10 -3 g cm -3 as the surface density value jPinzon fc Calvo-Mozol 
1200 ll) . The coefficient ipo was chosen in such a way that the maximum latitudinal velocity 
at middle latitudes is 20 m s _1 . The other values in eq. @ are: /3i=f.65 x fO 10 cm" 1 , 
/3 2 =2.2x fO 10 cm" 1 , e=2.0000001, r = {R Q - R min )/A.l 5, and T=3.47 x 10 10 cm. Here, R min 
is the minimum r coordinate value for t he integration rang e, and the free parameter Rb is the 
maximum depth of the return flow (see iM. Dikpatil If 9941 for more details). 



2.3 Magnetic buoyancy (MB) and the alpha effect 

As suggested in previou s simulations of magnetic flux tubes iS. D'Silval 1 1 9931 . |Y. Fanl Il993t 
IP. Caligarl Il995l Il998l) . a tube of strong magnetic field surrounded by a diffuse field can 
become buoyant unstable and arise to the surface where it is twisted by the coriolis force. It 
has been shown that a I0 5 G magnetic field becomes unstable and emerges to the surface in 
time intervals of about one month to form bipolar regions with the appropriate tilt angles. 
We introduce this effect in a simplified way by using a non-local regeneration of the toroidal 
magnetic field followed by an alpha effect which is concentrated close to the surface, so that: 

S 1 (r,9;B) = 

x 

X 

Where the parameters are r2=0.95i?Q, V3=Rq, and ^2=^3=0.025/?© . 



^B(r c ,e,t)[l + erf( r -^)] (10) 
[1 - erfC-^)] 

"3 

[1 + ( B^^ f]- 1 cos 9 sine . 
Bo 



2.4 Magnetic diffusivity 

The magnetic diffusivity is larger in regions of stronger magnet ic field. Therefore, we a ssume 
a larger diffusivity for the poloidal component, as suggested bv lChatteriee et all (l2004h : 

V P (r)=Vc+^[l + erf(2^^)] , (II) 

with r] c =2.2 x 10 9 and r] cz =0.2 x 10 11 cm 2 s _1 , and r c =O.77i0, and smaller diffusivity for the 
toroidal component: 

^M=^ + ^[l + ^/(^-p i )] + ^[l + er/(^£)] , (12) 

where jy cz i=4 x 10 10 cmV, r cl =0.72R Q , and r c2 =0.95i? Q . 
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3 Results 



We have solved the equations above by means of the ADI method in a two-dimensional mesh 
of 128x128 spatial divisions, with 0.55i?© < r < IRq, and < 6 < tt/2. In the model of Fig. 
1, the meridional flow is allowed to penetrate down to the base of the solar convection zone 
(r=0.675i? Q ). 
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Figure 1: Butterfly diagrams for (A) the radial field at the surface (r=R Q ) and (B) the toroidal 
magnetic field at the base of the solar convection zone (t^0.7Rq). The contours are equally 
spaced, with solid (dashed) lines for positive (negative) values. Time is in years and the latitude 
is in degrees. The maximum toroidal magnetic field at the base of the SCZ is located at latitudes 
between 75° and 90° with weak branches migrating toward the equator. With this result and based 
on magnetic flux tubes simulations, sunspots will appear mostly in regions close to the poles. 

3.1 Deep Meridional Flow 

If we allow the flow to go deeper (r=0.61i?e), a better distribution of the toroidal field lines 
is obta ined, however, a pola r bran ch also appears, in contrast to the observations and to the 
results iNandv fc Choudhuril J2002I) (Fig. 2). 
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Figure 2: The same as in Fig. 1, but for a deeper meridional flow. If most of the magnetic field 
is dragged below the tachocline, only a small portion emerges to the surface by the action of the 
magnetic buoyancy force. 
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3.2 Exploring an ellipsoidal shape of the tachocline 



Another way to test the deep meridional flow hypothesis is to change the form of the tachocline. 
In the helioseismic results of lCharbonneau et al.1 lll999f) . they have reported a possible prolate 
shape of the tachocline. Though the full properties of the tachocline still remain unknown, we 
may here modify the shape of this layer by replacing the constant term r c in the equations 
above for a term that depends on r and 9 r c (r,6): 



1 + ri v sin^ 



(13) 



where r cx and r cy are the semi-axis of an ellipse. If we first vary the term r c only in 
the differential rotation eqs. (6) and (10), in order to assure a substantial vertical shear, but 
maintaining it constant in eqs. (11) and (12), this configuration produces a prolate tachocline 
(r cx =0.67R.Q and r C y=0.73R,Q) that generates weak polar branches (with magnetic fields of 
the order of 10 4 G) and strong equatorward branches of toroidal field, in good agreement with 
the observations (see Fig. [3J|. In this case, the poloidal magnetic field, which is advected by 
meridional flows, goes deeper at the high latitudes within the more stable radiative layer. The 
toroidal field will diffuse more rapidly due to the longer values of the coefficient of diffusion 
in these latitudes. Advection then makes the rest of the process, draging this field to reach 
the tachocline and the convection zone, and arise to the surface. This scenario at which the 
flux goes deeper is more reliable because the mixture of elements occurs in a thinner layer 
within the radiative zone. If we now allow the r c term to vary according to eq. (13) also 
in the equations (11) and (12) and test this prolate configuration, we recover the results of 
the previous section (Fig. 2). A more detailed analysis will be presented elsewhere in a 
forthcoming paper. 




Figure 3: Butterfly diagrams for (A) the radial and (B) the toroidal magnetic fields assuming 
a prolate tachocline and a source term proportional to the toroidal magnetic field just above the 
tachocline. 



4 Conclusions 

We have recovered the iNandv fc Choudhuril i2002T) results and qualitatively reproduced the 
observed lati tudinal distribution of the mag netic fields in the Sun assuming a different source 
formulation jDikpati fc Charbonneaul FT9991 . The most important results here found may be 
summarized as follows: 

1. If the meridional flow is confined to the convective zone (r>0. 675Rq), the sunspots are 
mainly concentrated near the poles. 
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2. If the flow is allowed to penetrate deeper down to 0.61i?Q , the model gives a better result 
since a branch of maximum toroidal field is produced at low latitudes in agreement with 
the observations. However, in this case, another branch also appears near the poles which 
is inconsistent with the observed butterfly diagrams. 

3. When we incorporate a probably more realistic bipolytropic density profile in the tachocline 
and the radiative zone below the convective layer, the results above are maintained. 

4. Finally, if we consider an ellipsoidal tachocline instead of a spherical one, we find that 
when a prolate shape is assumed and the magnetic field is stored mostly above the 
tachocline, then the resulting butterfly diagram is in good agrement with the observa- 
tions. 

5. The results above suggest that, although the models with a deep meridional flow are, in 
general, able to produce results which are in good agreement with the observations, the 
global behavior of the circulation and the physical mechanism behind it needs further 
revision in order to determine where the magnetic field is really amplified and stor ed 
(see Guerrero, de Gouveia Dal Pino & Munoz 2005, in preparation). 
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